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We present a deep study of the average hard X-ray spectra of Seyfert galaxies. We aim to test the unified model of active galactic 
nuclei, and constrain differences and similarities between different classes of objects. We analyzed all public INTEGRAL IBIS/ISGRI 
data available on all the 165 Seyfert galaxies detected at z < 0.2. Our final sample consists of 44 Seyfert Is, 29 Seyfert 1.5s, 78 
Seyfert 2s, and 14 narrow-line Seyfert Is. For each subsample, we stacked all the images, and derived their average hard X-ray spectra 
in the 17-250 keV energy range. We performed a detailed spectral analysis using both a model-independent and a model-dependent 
approach. All classes of Seyfert galaxies show on average the same nuclear continuum, as foreseen by the zeroth order unified model, 
with a cutoff energy of Ec ^ 200 keV, and a photon index of T = 1.8. The average optical depth of the Comptonizing medium 
is consistent for the different classes (r a 0.8). Compton-thin Seyfert 2s show a reflection component stronger than Seyfert Is and 
Seyfert 1.5s. Most of this reflection is due to mildly obscured (10 23 cirT 2 < Nh < 10 24 cirT 2 ) Seyfert 2s, which have a significantly 
stronger reflection component (R = 2.2+ 4 ^) than Seyfert Is (if < 0.4), Seyfert 1.5s (R < 0.4) and lightly obscured (N a < 10 23 cnT 2 ) 
Seyfert 2s (R < 0.5). This cannot be explained easily by the unified model. The absorber/reflector in mildly obscured Seyfert 2s 
might cover a large fraction of the X-ray source, and contain clumps of Compton-thick material. The large reflection found in the 
spectrum of mildly obscured Seyfert 2s reduces the amount of Compton-thick objects needed to explain the peak of the cosmic X-ray 
background. Our results are consistent with the fraction of Compton-thick sources being ~ 10%. The spectra of Seyfert 2s with and 
without polarized broad lines do not show significant differences, the only difference between the two samples being the higher hard 
X-ray and bolometric luminosity of Seyfert 2s with polarized broad lines. The average hard X-ray spectrum of narrow-line Seyfert Is 
is steeper than those of Seyfert Is and Seyfert 1.5s, probably due to a lower energy of the cutoff. 

Key words. Galaxies: Seyferts - X-rays: galaxies - Galaxies: active - Galaxies: nuclei - X-rays: diffuse background 
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1. Introduction 

Active Galactic Nuclei (AGN) emit over the entire electromag- 
netic spectrum, and are commonly assumed to be powered by 
accretion onto a super massive black hole (Rees 1984). Seyfert 
galaxies host AGN, and are classified according to their emis- 
sion lines as Seyfert Is (Syls, showing broad and narrow emis- 
sion lines), and Seyfert 2s (Sy2s, showing only narrow lines). 
Many Seyfert galaxies exhibit optical spectra with properties 
in-between those of Syl and Sy2 galaxies (e.g., Osterbrock & 
Koski 1976). These objects have been classified as type 1.2, 1.5, 
1.8, or 1.9 Seyfert galaxies, depending on the details of their op- 
tical spectra. Studies of the optical spectra of Seyfert 2 galaxies 
using polarized light (Miller & Antonucci 1983, Antonucci & 
Miller 1985) showed that broad emission lines can also be de- 
tected in some of these objects, and led to the development of 
the so-called unified model (UM, Antonucci, 1993). According 
to this model, the same engine is at work in all kind of Seyfert 
galaxies, and the differences between Seyfert Is and Seyfert 2s 
can be ascribed solely to orientation effects and to anisotropic 
obscuration. In this model, the line of sight to the nucleus is 
(Seyfert 2s) or is not (Seyfert Is) obstructed by optically thick 
material, possibly distributed in a toroidal geometry. X-ray ob- 
servations have confirmed this idea, showing that most AGN un- 
absorbed in X-rays are of the optical Seyfert 1 type, and that 
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most AGN that are absorbed belong to the Seyfert 2 group (e.g., 
Awakietal. 1991). 

However, in the past few years observational evidence of sig- 
nificant differences between Seyfert Is and Seyfert 2s has been 
discovered. Seyfert Is with significant absorption have been 
found (e.g., Fiore et al. 2001, Cappi et al. 2006), along with 
Seyfert2s without X-ray absorption (e.g., Pappa et al. 2001, 
Panessa & Bassani 2002). Furthermore, spectropolarimetric sur- 
veys indicate that only ~ 30 - 50% of Seyfert 2s have polarized 
broad lines (PBLs), which might imply that not all of them har- 
bor hidden BLRs (e.g., Tran 2001, 2003, Gu & Huang 2002). 
In the soft X-rays the photon index distribution is found to be 
skewed towards low values of F for obscured objects (Brightman 
& Nandra 2010), while in the hard X-rays, Seyfert 2s have been 
found to have harder spectra than Seyfert Is (e.g., Zdziarski et 
al. 1995, Malizia et al. 2003, Deluit & Courvoisier 2003, Ajello 
et al. 2008a, Burlon et al. 201 1). 

Two subclasses of objects cannot be easily explained by the 
unified model: narrow-line Seyfert Is (NLSls) and low ioniza- 
tion nuclear emission line regions (LINERs). It has been sug- 
gested that NLSls are AGN in their early phase (Grupe et al. 
1999), characterized by relatively low black hole masses (e.g., 
Grupe & Mathur 2004) and very high accretion rates in terms 
of Eddington units (e.g., Grupe et al. 2010). LINERs (Heckman 
1980) are low luminosity AGN, and display multiwavelength pe- 
culiar characteristics. They could be the link between AGN and 
normal galaxies (Zhang et al. 2009, Rovilos et al. 2009), or their 
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SMBH might accrete differently from normal Seyfert galaxies 
(Ho 2008). 

The X-ray emission of Seyfert galaxies is thought to be pro- 
duced by the Comptonization of ultraviolet (UV) photons gener- 
ated in the innermost edge of the accretion disk by a population 
of hot electrons located in a coronal region sandwiching the disk 
(Haardt & Maraschi 1991, 1993). The X-ray spectrum of AGN 
can be normally well described by a power law model with a 
photon index between F ^ 1.8 (e.g., Dadina et al. 2008) and 
r =! 2 (e.g., Beckmann et al. 2009). Other common characteris- 
tics are an exponential cutoff at an energy depending on the tem- 
perature of the plasma T e and its Thomson opacity, and photo- 
electric absorption at low energies, caused by material along the 
line of sight. Two prominent features generated by the reflection 
of the continuum are also observed in the X-ray spectra of AGN: 
a neutral iron Ka line (at 6.4 keV in the local reference frame) 
and a reflection hump peaking at E ^ 30 keV (Magdziarz & 
Zdziarski 1995). At energies > 20keV, the photoelectric cross- 
section sharply declines, and the presence of matter along the 
line of sight does not play a significant role, unless the source 
is Compton-thick (CT; N H > 1.5 x 10 24 cnT 2 ). Therefore hard 
X-rays are ideal for probing the intrinsic emission of AGN and 
testing the UM, according to which different types of Seyfert 
galaxies should show on average the same characteristics. 

The INTErnational Gamma-Ray Astrophysics Laboratory 
(INTEGRAL; Winkler et al. 2003) is a hard X-ray /soft y-ray mis- 
sion, designed for imaging and spectroscopy with high angular 
and spectral resolution in the energy range from ~ 3 keV to 
10 MeV. Since its launch on October 17, 2002, the soft y-ray 
imager IBIS/ISGRI (15-1000 keV; Lebrun et al. 2003, Ubertini 
et al. 2003) on board INTEGRAL has detected more than 200 
AGN, among which about 90 were detected for the first time in 
the hard X-rays. IBIS/ISGRI uses the coded-aperture technique 
(Caroli et al. 1987) and has a large field of view of 29 ° square 
with a spatial resolution of 12 arcmin. 

We present a very deep study of the average hard X-ray spec- 
tra of Seyfert galaxies. Our main aim is to test the unified model 
by studying the average hard X-ray emission of different types 
of radio-quiet AGN detected by INTEGRAL IBIS/ISGRI, and 
to constrain the average spectral characteristics of each class. 
We also investigate whether Seyfert 2s showing (hereafter PBL 
Sy2s) and not showing (hereafter NPBL Sy2s) PBLs are intrin- 
sically similar objects or not. In Sect. [2] we present the sample 
of AGN; in Sect. [3] we describe how the data analysis was per- 
formed and how the hard X-ray spectra were obtained; in Sect.lU 
we present the results obtained from a model-independent analy- 
sis of the spectra, and in Sect. [5] we discuss the model-dependent 
spectral analysis. We study the origin of a strong reflection com- 
ponent found in the spectrum of Seyfert 2s in Sect. [6] and in 
Sect. 13 we discuss our results and their implications for the UM 
and the cosmic X-ray background (CXB). In Sect.|8] we present 
our conclusions and summarize our findings. 



2. The sample 

The sample consists of all the 205 Seyfert galaxies detected (sig- 
nificance > 5cr) in the hard X-rays by INTEGRAL IBIS/ISGRI 
during its first eight years of operations. We considered as their 
optical classifications the ones reported in the latest version of 
the Veron Catalog of Quasars and AGN (Veron-Cetty & Veron, 
2010). 

The sample consists of 38 Seyfert Is, 17 Seyfert 1.2s, 34 
Seyfert 1.5s, 3 Seyfert 1.8s, 15 Seyfert 1.9s, 75 Seyfert2s, 8 



LINERs, and 15 NLSls. We excluded from these subsam- 
ples the LINERs, because the physical characteristics of their 
accretion flow might differ from those of the others. We 
grouped Seyfert 1.2s with Seyfert Is, and both Seyfert 1.8s and 
Seyfert 1.9s with Seyfert 2s because of their similarity. 

2.1. The role of absorption 

To model the effect of absorption on the X-ray spectra of AGN, 
one has to take into account both photoelectric absorption and 
Compton scattering. The photoelectric cross-section cr ph has a 
strong dependence on the energy, while Compton scattering in 
the hard X-rays depends on the Thomson cross-section o"t, and 
is constant with energy up to the Klein-Nishina decline. The cu- 
mulative effect of the two cross-sections is given by 



(1) 



For column densities of the order of cr T ' = 1.5 x 10 24 cm" 2 , 
only about 30% of the flux is able to escape at ~ 20keV (see 
Fig. [I). About 60% of the original emission is unabsorbed for 
A^h = 7 x 10 23 cm -2 , and we used this value as a threshold be- 
tween Compton-thin and Compton-thick sources. Thus, we in- 
cluded in our final sample of Seyfert 2s only the objects with a 
value of column density of A^h < 7 x 10 23 cirT 2 , and placed 
the 10 Seyfert 2s with N H > 7 x 10 23 cnr 2 in the CT Seyfert 2s 
sample. We note, however, that none of the sources in our sam- 
ple has a A^h in the (7-10)xl0 23 crrT 2 range, so in the follow- 
ing we refer to 10 24 crrT 2 as the threshold between Compton- 
thin and Compton-thick objects. In the following, we use solar 
abundances and the photoelectric cross section of Morrison and 
McCammon(1983). 
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Fig. 1. Fraction of escaping flux in the 10 energy bins used. 
Squares represent A^h = 10 23 crrT 2 , stars A^h = 7xl0 23 cirT 2 , di- 
amonds Nu = 1 .5 xl0 24 cirT 2 , and triangles Af H =5xl0 24 cm~ 2 . 



2.2. Excluded sources 

To avoid contamination and problems related to the way we ob- 
tain the average spectra (see Sect. [3), we excluded from the sam- 
ple 10 sources which are too close (< 0.4°) to each other or 
to other bright sources. We excluded 1 1 sources with redshift 
z > 0.2 to avoid spurious effects due to redshifted spectral fea- 
tures (i.e. the high-energy cutoff and the reflection hump). We 
removed the 5 sources detected with a significance >50<r, to 
prevent them from dominating the final spectra. 

Four sources with either an uncertain classification or pe- 
culiar absorbers were also excluded from the analysis. Was 49 
is an interacting system of galaxies containing two Seyfert nu- 
clei (Moran et al. 1992). The AGN IGR J16426+6536 shows 
characteristics typical of both NLS1 (Masetti et al. 2009) 
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Table 1. Sources detected by INTEGRAL IBIS/ISGRI but ex- 
cluded from the analysis, and reason for their exclusion. More 
details can be found in Sect. 12.21 



Source 


Type 


Reason for the 


4TT H44-60 


Sy 1.5 


c 




CT Sy 2 


B 


FSO 1 38-1 


Sy 2 


c 


FSO 383-18 

1 J v.f V / JUJ x o 


Sy 2 


c 


GRS 1734-292 


Sy 1 


B 


tc 4329A 


Sv 1 
oy i 


B 


TOR T00465-4005 


Sy 2 




TGR T05270-6631 


Sy 1 




TGR T061 17-6625 

XVJXV JUU1 X 1 VJVJZ..J 


Sy 1.5 




TGR T09523-6231 

1 V J I V JU7JZJ UiJl 


Sy 2 




TGR T1 01 47-6354 


Sy 1 




TGR T1213 1+0700 

1 V I XV J 1 — 1 ^' 1 ~ V/ / \J\J 


Sy 1.5 


c 


TGR T1 3038+5348 

XVJXV J 1 JUJOT J JT^U 


Sy 1 


n 


TGR T1 45 15-5542 


Sy 2 


c 


TGR T 16024-6 107 


Sy 2 


c 


TGR T 16056-61 10 

XVJ XV J 1UWJU Ul lu 


Sy 1.5 


c 


TGR T1 6426+6536 

XVJXV J X VJ^ ^ VJ 1 VJ .J .J W 


NLS1 


U 


TGR T1 7476-2251 

XVJXV J 1 1 1 / VJ Z.Z.J J 


Sy 1 


c 


IGR J17488-3253 


Sy 1 


c 


IGR J 18249-3243 


Sy 1 


R 


IGR J2 1272+4241 


Sy2 


R 


NGC 1365 


Sy2 


U 


NGC4151 


Sy 1.5 


B 


NGC 4388 


Sy2 


B 


NGC 4939 


Sy2 


U 


NGC 6221 


Sy2 


C 


PKS 0637-752 


Sy 1.5 


R 


QSOB 182 1+643 


Sy 1 


R 


SWIFT J0216.3+5128 


Sy2 


R 


Was 49a + Was 49b 


Sy 1 + Sy 2 


U 



Notes. C: close to other detected sources (< 0.4°); U: classification 
unclear or with particular characteristics; R. z> 0.2; B: bright 
(detection significance > 50cr in the 17-80 keV band). 
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Fig. 2. Hydrogen column density distributions of Seyfert Is, 
Seyfert 1.5s, Seyfert 2s and narrow-line Seyfert Is. 



and Seyfert 1.5 (Butler et al. 2009) galaxies. The Seyfert 1.8 
NGC 1365 has been found to undergo transitions from a "trans- 
mission dominated" to a "reprocessing dominated" state and 
back in a few weeks time, which have been interpreted as due 
to eclipses of CT material along the line of sight (Risaliti et al. 
2005). A similar state transition has been observed in NGC 4939 
(Guainazzi et al. 2005), a Seyfert 2 that was previously classi- 
fied as Compton-thick by BeppoSAX observations (Maiolino et 
al. 1998). The excluded sources are reported in TableQ] 

2.3. The final sample 

The final sample used for the analysis consists of 165 AGN, 
of which 44 are Seyfert Is, 29 Seyfert 1.5s, 68 Compton-thin 
Seyfert2s, 10 CT Seyfert 2s, and 14 NLSls. Spectropolarimetric 
data are available for only 19 of the 68 Compton-thin Seyfert 2s 
of our sample. Among these, 1 1 show and 8 do not show PBLs. 
The sources and their characteristics are reported in AppendixlAl 

Figure|2] shows the distribution of the column density for 
the different classes of AGN. The average values of the column 
densities of the different classes are A^ yl = 5.4 x 10 21 cm , 
N S H y15 = 3.7 x 10 22 cirT 2 , N S R y2 = 1.5 X 10 23 cirT 2 , A^ LSI = 
1.9 x 10 21 cirT 2 , and A^ T = 1.5 x 10 24 cirT 2 , for Seyfert Is, 
Seyfert 1.5s, Seyfert 2s, NLSls, and CT Seyfert 2s, respectively. 
These values agree with the paradigm that Seyfert 2s are more 
obscured than Seyfert Is and Seyfert 1.5s. Considering the aver- 
age column densities of the different samples, the percentage of 
flux absorbed in the first band (17-22 keV) by photoelectric pro- 
cesses is about 0.01%, 0.5%, and 2% for Seyfert Is, Seyfert 1.5s, 
and Seyfert 2s, respectively. Compton processes have a stronger 
influence at these energies, and the fraction of scattered flux is of 
0.3%, 2.4% and 9.5%, for Seyfert Is, Seyfert 1.5s and Seyfert 2s, 
respectively. In the case of CT Seyfert 2s, the influence of the ob- 
scuring material is stronger and 18% of the flux is photoelec- 
trically absorbed, whereas when we also consider Compton scat- 
tering, only 30% of the flux escapes. A Kolmogorov-Smirnov 
(KS) gives a probability of =s 67% that the column density dis- 
tributions of Seyfert 1 s and Seyfert 1 .5s are drawn from the s ame 
parent population. The probability is much lower when compar- 
ing NLSls to Seyfert Is and Seyfert 1.5s (=* 10%). 

In Fig. [3] we show the redshift distribution of our final sam- 
ple. The average value of the redshift is z=0.03, and the average 
detection significance is 12. Icr. 

In Fig. |H we show the luminosity distributions in the 17- 
80 keV band of our samples of Seyfert galaxies. The luminosi- 
ties were calculated assuming a standard ACDM cosmological 
model, Crab-like spectra (e.g., Jourdain & Roques 2009), and 
Ho = 70kms~ 1 Mpc -1 . The average values of the luminos- 
ity are L Sy i = 5.0 x lO^ergs -1 , L Sy i.5 = 3.7 x lO^ergs" 1 , 
L Sy2 = 2.1 x 10 43 ergs _1 , L CT s y 2 = 8.3 x 10 42 ergs _1 , and 
£nlsi = 1.7 xlO 43 ergs" 1 , for Syls, Syl.5s, Sy2s,CTSy2s, and 
NLSls, respectively. The luminosity of CT Sy2 was corrected 
for Compton scattering. A KS test gives a probability of =2 53% 
that the luminosity distributions of Seyfert Is and Seyfert 1.5s 
are statistically compatible. The probability is much lower when 
comparing Seyfert Is to Seyfert2s (= 2%), and Seyfert Is to 
NLSls (=10%). 

3. INTEGRAL IBIS/ISGRI data analysis 

We used all the public data obtained by INTEGRAL IBIS/ISGRI 
as of May 2010, for a total of more than 50 thousands point- 
ings or "science windows" (SCWs). The typical exposure time 
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Fig. 3. Redshift distribution of our final sample of INTEGRAL 
IBIS/ISGRI detected Seyfert galaxies. 
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Fig. 4. Luminosity distribution of Seyfert Is, Seyfert 1.5s, 
Compton-thin Seyfert 2s, NLS Is, andCTSy2s in the 17-80 keV 
energy band. 



of these pointings is of (1 - 3) X 10 3 s. Several thousand point- 
ings, including at least one source of the sample in its field of 
view and with an effective exposure longer than 120 s were se- 
lected, spanning times between December 30, 2002 (revolution 
26) and April 7, 2009 (revolution 791). 

The ISGRI data were reduced using the INTEGRAL Offline 
Scientific Analysis softwar43 version 9.0, publicly released by 
the ISDC Data Centre for Astrophysics (Courvoisier et al. 2003). 

Individual sky images for each pointing were produced in 
a broad energy band (17-250 keV), and divided into ten bins 
as follows: 17-22, 22-30, 30^0, 40-51, 51-63, 63-71, 71-80, 
80-105, 105-150 , and 150-250 keV. 

To extract the average spectra of the different AGN samples, 
we followed the procedure adopted by Walter & Cabral (2009). 
We created 500x500-pixels mosaic images modifying the coor- 
dinate system of each individual image, setting the coordinates 
of each source of the sample to an arbitrary fixed position (or=0, 
6=0). The geometry of the image was also modified to have a 
consistent PSF independently of the position of the source in the 



Table 2. Number of sources (1), effective exposures (2), detec- 
tion significances (3), and count rates (4) of the different sam- 
ples. All the values refer to the 17-80 keV band. 







\A) 




w 


Sample 


Srcs 


Exp. 


Det. Significance 


Ct rate 






[Ms] 


H 


[ct/s] 


Seyfert 1 


44 


37.6 


91.1 


0.396 ± 0.004 


Seyfert 1.5 


29 


13.2 


59.1 


0.429 ± 0.007 


Seyfert 2 


68 


44.0 


102.8 


0.390 ± 0.004 


CT Sy2 


10 


4.7 


32.9 


0.405 ± 0.012 


NLS1 


14 


10.1 


44.3 


0.377 ± 0.008 


MOB Sy2 


27 


13.8 


59.6 


0.430 ± 0.007 


LOB Sy2 


34 


23.4 


78.1 


0.424 ± 0.006 


PBL Sy2 


11 


4.2 


99.0 


1.301 ±0.013 


NPBL Sy2 


8 


2.9 


31.2 


0.479 ± 0.015 



field of view (FOV). These mosaic images provide a stack of all 
the selected IBIS/ISGRI data for each considered sample. 

To minimize the systematics in the mosaic images, we ex- 
cluded 1455 individual sky images (many obtained before revo- 
lution 38, when the IBIS bottom anti-coincidence was reconfig- 
ured) that had a background fluctuation rms larger than 1 Act in 
the significance image, and 779 images with a minimum signifi- 
cance smaller than -5.5o\ In addition to these, we excluded 627 
images taken when the ISGRI was in staring mode. A total of 
56,61 1 images were finally included in the processing (some of 
them many times, when including several sources). All known 
ISGRI sources with significance above 5<x were used for image 
cleaning. When known accurately, the position of these sources 
was fixed to catalogue values taken from the INTEGRAL general 
reference catalog (Ebisawa et al. 2003 

The mosaic images were built with a tangential projection 
using a factor of two oversampling when compared to the indi- 
vidual input sky images, this results in a pixel size of 2.4 arcmin 
at the center of the mosaic. The photometric integrity and ac- 
curate astrometry were obtained by calculating the intersection 
between input and output pixels, and weighting the count rates 
according to the overlapping area. 

The average signal extracted in the large 17-80 keV band 
from each ISGRI mosaic, and both the exposures and the num- 
ber of sources used are reported in Table |2] The detection sig- 
nificances in the 17-80 keV band range between 102. 8cr and 
32.9<r, for the largest and smallest sample (i.e. Seyfert 2s and CT 
Sy2s), respectively. The effective exposure obtained at the cen- 
ter of the mosaics are between 4.7 Ms (for CT Sy2s) and 44.0 Ms 
(for Seyfert 2s). 

The ten-bin spectra of the final samples were extracted from 
the images in the 10 narrow bands using mosaic_spec. In all 
spectra, the highest energy bin (i.e. 150-250 keV) has a low sig- 
nificance (~ 2.5 - 3cr), but we use it for the sake of complete- 
ness. We used the latest detector response matrix (RMF), and 
calculated the ancillary response matrices (ARFs) of each sam- 
ple with a weighted average of the nine available ARFs, based 
on the number of SCWs within the validity time of a particular 
ARF 



1 http://www.isdc.unige.ch/ 



2 http : //www . isdc . unige . ch/integral/ science/catalogue 
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4. Model-independent spectral analysis 

A model-dependent spectral analysis has the drawback of pos- 
sible parameter degeneracy, which might not allow us to con- 
strain several parameters at the same time. An alternative 
method to characterize differences and similarities between dif- 
ferent classes of Seyfert galaxies, independently of their average 
flux, is by means of a model-independent approach. This was 
achieved by normalizing the flux of the two different spectra in 
the first bin (17-22 keV), and then calculating their ratio. 



4.1. Seyfert galaxies 

We have no information about the value of the hydrogen column 
density of seven Seyfert 2s. To discard the possibility that these 
objects are CT, which might influence significantly the result, we 
calculated the ratio between the spectra of the Seyfert 2 sample 
including and excluding these seven objects. This ratio is fully 
consistent with 1 . 

In Fig. [5] we show the normalized spectra of Syls and Syl.5s 
(in the upper panels) and their ratios (in the lower panel). 
The average spectrum of Seyfert Is is consistent with that of 
Seyfert 1.5s along the whole spectrum. The ratio of the average 
spectrum of Seyfert 1 galaxies to the normalized one of NLSls 
diverges from the unity from ~ 40keV (Pig.[6j, with the spec- 
trum of NLSls being steeper than that of Syls. Since the spec- 
tra of Seyfert Is and Seyfert 1.5s are very similar, we merged 
them and compared the resulting spectrum to that of Seyfert 2s. 
We show in Fig. [7] the ratio obtained after normalization. The 
spectrum of Seyfert 2s displays a clear excess over those of the 
Seyfert Is and Seyfert 1.5s of 48 + 8% in the 22-63 keV band. 
This bump peaks around 40keV, reaching 20%, and might in- 
dicate that Seyfert 2s have a stronger reflection component than 
Seyfert Is and Seyfert 1.5s. At energies higher than 60keV, the 
ratio is consistent with 1, indicating similar cutoff energies. The 
ratio of the spectra of Compton-thin to Compton-thick Seyfert 2s 
is consistent with one up to 150 keV (Fig. [8]). Including in our 
sample of Sy2s the two objects excluded because of their com- 
plex absorption does not affect the results, and the ratio between 
the spectra with and without these sources is consistent with one. 



Seyfert 1 vs Seyfert 1.5 




00 ISO 200 750 



Fig. 5. Normalized spectra (upper panel) and ratio (lower panel) 
of Seyfert Is to Seyfert 1.5s. 




Seyfert 1 vs NLS1 



100 150 200 250 



Fig. 6. Normalized spectra (upper panel) and ratio (lower panel) 
of Seyfert Is to NLSls. 




Energy (keV) 



Fig. 7. Ratio between the normalized spectra of Compton-thin 
Seyfert 2s and that of both Seyfert Is and Seyfert 1.5s. 



4.2. Lightly obscured and mildly obscured Seyfert 2s 

The Seyfert 2s contained in our sample appear to have char- 
acteristics that differ significantly from those of Seyfert 1 and 
Seyfert 1.5 galaxies. To test whether absorption might still play 
a role in the hard X-ray spectra of Compton-thin Seyfert 2 galax- 
ies, we divided the sample into two subsamples, one includ- 
ing the objects with 10 23 cnf 2 < Nh < 10 24 cm~ 2 (henceforth 
mildly obscured Seyfert 2s, MOB Sy2s), and the other those 
with jVh < 10 23 cirT 2 (henceforth lightly obscured Seyfert 2s, 
LOB Sy2s). The MOB and LOB Sy2 samples contain 27 and 35 
objects, respectively. In Table|2] we report the detection signifi- 
cances, exposures, and fluxes obtained for these samples. We did 
not consider the sources with an unknown value of A^h- 

In Fig. |9l we show the ratio obtained by comparing the spec- 
trum of MOB Sy2s to that of LOB Sy2s. The figure highlights 
significant differences between the two spectra, with the emis- 
sion of MOB Sy2s being harder up to ~ 60 keV and then softer 
(up to ~ 150 keV) than that of LOB Sy2s. This clearly shows 
that the bump seen by comparing the spectrum of Seyfert 2s to 
those of Seyfert Is and Seyfert 1.5s (Fig. [7) comes prevalently 
from the contribution of MOB Sy2s (see Fig.lTOb. 

The spectrum of LOB Sy2s appears to be very similar to that 
of Seyfert Is (Fig.lTOt. although an excess of ^ 20% is evident in 
the 30-50 keV band. This excess might be related to a stronger 
reflection component in LOB Seyfert 2s than in Seyfert Is. MOB 
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Sy2s present instead a spectrum similar to that of CT Sy2s, and 
their ratio is consistent with 1 (Fig.fTOll. 

To test the influence of bright MOB Sy2s on their average 
spectrum, we compared the spectra obtained with and without 
the five brightest sources (> 20<x). The ratio obtained after nor- 
malizing the two spectra is fully consistent with one, from which 
we can conclude that the average spectrum of MOB Sy2s is not 
influenced by peculiar bright sources. 



CT Sy2 vs Seyfert 2 




150 200 250 



Fig. 8. Normalized spectra (upper panel) and ratio (lower panel) 
of Compton-thin Seyfert 2s to CT Sy2s. 
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MOB Sy2 vs CT Sy2 
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Fig. 10. Ratio between the normalized spectra of LOB Seyfert 2s 
and Seyfert Is (top panel), between the spectra of MOB 
Seyfert 2s and Seyfert 1 s (center panel), and between the spectra 
of MOB Seyfert 2s and CT Seyfert 2s (bottom panel). 



NPBL Sy2 vs PBL Sy2 
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Fig. 9. Ratio of the normalized spectra of MOB to LOB Seyfert 2 Fi §- ll - Ratio between the normalized spectra of PBL and NPBL 

galaxies Sy2s (top panel), between Seyfert Is and PBL (center panel), 

and NPBL Sy2s (bottom panel). 



4.3. Seyfert 2s with and without PBLs 

To test whether PBL and NPBL Sy2s have different characteris- 
tics, we extracted and analyzed their average hard X-ray spectra. 
In Table|2l we report the detection significances, exposures, and 
fluxes obtained for the two samples. 

Comparing the spectrum of PBL Sy2s to the one of NPBL 
Sy2s (Fig.fTTTi. it is evident that the ratio is consistent with 1. 

When compared to the spectrum of Seyfert Is, the spectra 
of PBL and NPBL Sy2s have similar characteristics. As for the 
whole Seyfert 2 sample, both PBL and NPBL Sy2s display evi- 
dence of a greater reflection component. This is probably due to 
the contribution of the MOB Sy2s present in the samples. 



4.4. Comparison with a simulated power law spectrum 

In Fig. Q2] we show the ratios between the normalized spectra 
of Seyfert galaxies and a simulated power law spectrum with a 
photon index of T — 1.9. Using XSPEC (version 12.5.0, Arnaud 
1996), we simulated 20 power-law spectra, using the fake it 
command. We used the latest RMF and calculated the ARF from 
a weighted average of the nine available ARFs, with the weights 
proportional to the validity period duration of each ARF. The 
power law spectrum used is the average of the 20 simulated spec- 
tra, and the errors are their standard deviations. 

Seyfert Is and Seyfert 1.5s show a ratio close to unity below 
100 keV, and a small excess around 30 keV, which is possibly 
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a signature of reflection. The excess in the 22-40 keV band is 
8 ± 3% and 9 + 3% for Seyfert Is and Seyfert 1.5s, respectively. 

The spectrum of Seyfert 2s displays a stronger excess (27 + 
3%) in the same band, and is harder than the power law up to 80 
keV. Above 150keV, the spectra of Seyfert Is, Seyfert 1.5s and 
Seyfert 2s become softer than the power law; this effect might be 
related to the presence of a high-energy cutoff. 

LOB Sy2s show an excess of 17 + 3% over the simulated 
power law in the 22-40 keV band, and do not present any signif- 
icant softening above 100 keV. MOB Sy2s show the most signif- 
icant differences, being harder below 60 ke V and softer above 
80 keV. The excess is of 87 + 9 % in the 22-63 keV band, and 
of 34 + 5 % in the 20^10 keV band. A similar excess (78 + 17% 
and 36 + 9%) is found by comparing the average spectrum of 
CT Seyfert 2s to the simulated power law in the same bands. 

To estimate the amount of curvature in the spectrum, we in- 
troduce the parameter p R given by 



fc ; Seyfert 1 



p R =1.09 ± 0.07 



PR 



^30 

Rgo 



(2) 



where R30 and R^q are the ratios of the normalized spectra to the 
power law in the 30-40 keV and 80-105 keV band, respectively. 
Seyfert Is (p R = 1.09 + 0.07) have a value of p R consistent with 
those of Seyfert 1.5s (p R = 1.06 + 0.11) and LOB Sy2s (p R = 
1.14 + 0.03). The value of this parameter is significantly larger 
for MOB Sy2s (p R = 1.63 + 0.14), while it is p R = 1.27 + 0.24 
for CT Sy2s. 



5. Model-dependent spectral analysis 

We analyzed the average IBIS/ISGRI spectra using XSPEC. 
Throughout the paper, we used 3<x errors. The probabilities that 
additive parameters improve a fit are calculated with the F-test. 

We tested simple models such as a power law, and a power 
law with an exponential cutoff at high energies. To account for 
the presence of reflection, we used the pexrav model (Magdziarz 
& Zdziarski 1995). Pexrav calculates the spectrum produced by 
an X-ray source reflected by an optically thick neutral slab. In 
this model, the parameter R measures the reflection component. 
If the source is isotropic, R can be linked to the solid angle Q 
as/? 



In 



The value of R depends on the inclination angle i 
between the perpendicular to the accretion disk and the line of 
sight. The smaller the inclination angle, the larger the resulting 
reflection component. The quality of the data does not allow us 
to constrain both R and i at the same time, so we fixed the in- 
clination angles to i = 45° for all the classes. We used the same 
value of i for all the classes to characterize the reflection with 
a single parameter. We note that by doing this R no longer rep- 
resents Q/2n, and becomes only an indicator of the reflection 
amplitude. 



5.1. Seyfert galaxies 

Fitting the spectra of Seyfert Is, Seyfert 1.5s, and Seyfert 2s with 
a simple power law yields a reduced chi-squared of xl - 1 -4, 
xl — 1.4, and^ 3.9, respectively. The null-hypothesis prob- 
ability is 20%, 18%, and < 0.1% for Syls, Syl.5s, and Sy2s, 
respectively. 

Applying a cutoff power law improves the fit of both 
Seyfert Is and Seyfert 2s at a confidence level of 84% and 99.6%, 
respectively. Adding a reflection component to the power law 
model also significantly improves the fit for both Seyfert Is and 
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Fig. 12. Ratio between the normalized spectra of Seyfert Is, 
Seyfert 1.5s, Compton-thin Seyfert 2, LOB Seyfert 2s, MOB 
Seyfert 2s, and Compton-thick Seyfert 2s, and a simulated power 
law spectrum (F = 1.9). The parameter p R is defined in Eq.[2] 



Seyfert 2s. The fit to the spectrum of Seyfert 1.5s does not im- 
prove significantly by adding any of these features to the baseline 
power law model. Fitting with different models the spectrum of 
Seyfert 1.5s, one notes that most of the chi-squared comes from a 
single bin (71-80keV, see AppendixIBl. On the basis of this, and 
that in Sect. |4] we have shown that the spectra of Seyfert Is and 
Seyfert 1.5s are very similar, we use the same models to char- 
acterize their spectra. The results of the spectral analysis are re- 
ported in Table|3] 

Both a high-energy cutoff and a reflection component from 
neutral matter are thought to play an important role in the hard 
X-ray emission of Seyfert galaxies, thus in the following we 
discuss only the results obtained using pexrav. Applying this 
model, we obtained a lower limit to the high-energy cutoff of 
Ec Z 190 keV, and a reflection normalization of R < 0.4 for both 
Seyfert Is and Seyfert 1.5s. The consistency of these parameters 
agrees with that found in Sect.|4] Using the same model for the 
average spectrum of Seyfert 2s, we obtained Ec — 154^ 5 keV 
and R < 2.6, which confirms the greater curvature observed. 
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Owing to its low statistics, a simple power law model pro- 
vides a good fit (xl = 1) to the stacked spectrum of the ten 
CT Seyfert 2 galaxies of our sample. Although statistically not 
required, a cutoff and a reflection component reduce the chi- 
squaredof A^ 2 =3.8 and Ax 2 = 5 . 1 , respectively. Using pexrav, 
we could not constrain Ec, while we obtained an upper limit to 
the reflection normalization of R < 1 1.2. 

A power-law model cannot represent well the hard X-ray 
spectrum of NLSls. Adding a cutoff or a reflection component 
significantly improves the fit, although it is not possible to con- 
strain the two parameters at the same time. We obtained a photon 
index of T = 2.3*2'i, a lower limit to the energy of the cutoff of 
Ec > 53 keV, and to the reflection parameter of R > 0. 1 . 

The average photon indices obtained are consistent for the 
different subsamples of Seyfert galaxies, and have values of T ^ 
1 .8, in agreement with previous studies (e.g., Dadina et al. 2008). 
The spectra of Seyfert Is, Seyfert 1.5s, Seyfert 2s, CT Seyfert 2s, 
and NLS 1 s are shown in AppendixlBl 

5.2. MOB and LOB Seyfert2s 

Using pexrav to fit the spectra of MOB and LOB Sy2s, we con- 
firmed the different curvatures of the spectra of the two sub- 
samples, obtaining values of the reflection parameter of R < 
0.5 and R = 2.2+^, for LOB and MOB Sy2s, respectively. 
The photon indices (F ~ 1.8) and the energies of the cutoff 
( £LOB _ 425+2« keV and £ mob _ 287+JfkeV) obtained by 

the fit are consistent to within 3<x . The results are reported in 
Table|4] while the spectra of MOB and LOB Sy2s are shown in 
AppendixlBl 

The spectral analysis performed using different models con- 
firms that the average hard X-ray spectrum of LOB Sy2s is more 
similar to those of Seyfert Is and Seyfert 1.5s than to that of 
MOB Sy2s. 

5.3. Seyfert 2s with and without PBL 

Consistently with what was obtained for the whole sample of 
Seyfert 2s, the spectra of PBL and NPBL Sy2s cannot be repro- 
duced using solely a power law, but need to be fitted using a 
more complex model. Adding a high-energy cutoff or a reflec- 
tion component to the power law significantly improves the fit 
for both spectra. All the models used give values of the param- 
eters of the two spectra that are consistent to within 3<x, con- 
firming the similarities found by the model-independent analy- 
sis. The results of the spectral analysis are reported in Table|5] 

6. Reflection in MOB Sy2s 

We have shown that the spectrum of MOB Sy2s shows a stronger 
reflection component than those of LOB Sy2s, Seyfert Is, and 
Seyfert 1.5s, and we now investigate the possible causes of this 
effect. 

The spectrum of lightly obscured objects (LOB Sy2s, Syls, 
and Syl.5s) can be expressed by 

F(E) = AE- T e- E/Ec + R(E), (3) 

where A is the normalization and R(E) the reflection component, 
which was accounted for using pexrav. The observed 20-60 ke V 
bump might be due to a distant (i.e. unabsorbed) reflection com- 
ponent, which becomes more important for MOB Sy2s due to 
the larger fraction of Compton scattered continuum. This reflec- 
tor has been associated with the inner wall of the torus, or the 



Table 3. Results obtained from the spectral analysis of the aver- 
age hard X-ray spectra of Seyfert Is, Seyfert 1.5s, Seyfert2s, CT 
Sy2s, and NLSls. The inclination angle i was fixed to 45°. 



Model 




r 


E c 


R 


X 2 /DO¥ 


FTEST 








[keV] 








Seyfert Is 


Power law 




1 Q7+0.05 
Ly -0.05 






11/8 




Cut-off power 


law 


1 s+°i 

L -°-0.2 


9CH+JVC 
• LyJ -150 




8.2/7 


0.17 


'Pexrav 








2 +a2 

"•^-0.2 


9.8/7 


0.38 


Pexrav 






J -155 


l +03 

"•4_0J 


8.3/6 


0.43 


Seyfert 1.5s 


Power law 




i Q7+0.09 
1,7 -0.08 






11.4/8 




Cut-off power 


law 


1 9 +01 


500 +31 ° 

Ju -250 




11.3/7 


0.81 


'Pexrav 




1 Q7+0.05 
Ly -0.05 




2 +a3 
" -0.2 


11.3/7 


0.81 


Pexrav 






T-59+iVC 
-145 


l +a3 

"•4_0J 


11.3/6 


0.97 


Seyfert 2s 


Power law 




1 QQ+0.03 

1 ' yJ -0.03 






31/8 




Cut-off power 


law 


1 4+0.2 






8.9/7 


<0.01 


'Pexrav 




1 Q7+0.05 
7 -0.05 




2 +3 9 


9.3/7 


< 0.01 


Pexrav 




1 6 +0/t 
'■•"-OA 




4 +22 
u - -0.4 


7.4/6 


0.01 


CT Seyfert 2s 


Power law 




1 9 +01 






8/8 




Cut-off power 


law 


1 5 +a2 


lOO^f 




4.2/7 


0.04 


'Pexrav 




1 9 +0 1 




1 4 +96 


2.7/7 


<0.01 


Pexrav 




2 +al 


NC 


1 5 +9J 

-1.5 


2.7/6 


0.04 


Narrow Line 


Seyfert Is 










Power law 




9 9 Q +0.08 
•"•"-0.08 






24.6/8 




Cut-off power 


law 


1 7 +0 - 2 

1 ■ ' -0.6 


7Q+43 
'"-36 




11.6/7 


0.03 


'Pexrav 




9 90 +0.08 
2 - ZS -0.08 




4 T.+NC 
^ -3.0 


9.5/7 


0.01 


Pexrav 




9 O+0.1 

^■ J -0.6 


3,0+jvc 

-257 


4 2 +NC 
^- -4.1 


9.5/6 


0.06 



Notes. NC: the parameter or the 3cr error were not constrained. 
'Pexrav: the high-energy cut off was not included, being fixed to 
its upper limit E c = 10 6 keV. 



outer part of the disk. The spectrum of MOB objects can be rep- 
resented by 

F(E) = M(E, Nn)AE- r e- E/Ec + R(E). (4) 

In the following, we assume that F = 1.8, Ec = 300 keV. We 
first simulated and compared several spectra with the same con- 
tinuum and reflection component, but with different values of 
/Vh, to see whether absorption might explain the differences ob- 
served. The value of R associated with the torus remains un- 
known (e.g., Gilli et al. 2007 consider a value of Rj = 0.37). 
To obtain an upper limit to the effect, we assumed R — 1 (for 
both LOB and MOB objects), and column densities of /Vh = 
7 x 10 23 cirT 2 and 7Y H = for MOB and LOB objects, respec- 
tively. We tested two different values of the inclination angle for 
MOB Sy2s (i = 30° and 2 = 60°), while we set i = 60° for 
LOB objects. The parameter M(E,Nh) includes both Compton 
scattering (through the cabs model) and photoelectric absorp- 
tion (implemented as wabs), as given in EqQ] In Fig. [13] {left 
panel), we show the ratios obtained comparing the simulated 
spectra. From the figure, it is evident that a column density of 
/Vh = 7 x 10 23 ctrT 2 is largely insufficient to explain the ob- 
served differences. To provide an excess closer to what is ob- 
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Fig. 13. Ratios of simulated spectra of MOB (Eq.|4ji to LOB (Eq.0 objects. We assumed r = 1.8, Ec — 300 keV for the continuum 
of all spectra. We tested different values of i and A^h (A^h = 7 x 10 23 crrT 2 and A^h = 5 x 10 24 cm -2 ) for the spectrum of MOB Sy2s, 
assuming R = 1 for all the spectra, and A^h = for LOB Sy2s (left panel). We tested different values of the reflection parameter 
(R = 0.4 and R = 0.2 for LOB objects, fixing R = 1 f or MOB Sy2s), and of i (only for MOB Sy2s), setting N H = 3 x 10 23 cm -2 for 
MOB Sy2s and A^h — for LOB Sy2s (right panel). 



Table 4. Results obtained from the spectral analysis of the aver- 
age hard X-ray spectra of MOB and LOB Seyfert 2s. The incli- 
nation angle i was fixed to 45°. 



Model 


r 


E c 


R 


* 2 /DOF 


FTEST 






[keV] 








MOB Seyfert 2s 


Power law 


1 Q4+0.05 
• -0.05 






45.2/8 




Cut-off power law 


1 +03 

^ -0.3 


47 +20 




7.6/7 


< 0.01 


'Pexrav 






7 +7 
-2 


9.8/7 


< 0.01 


Pexrav 


1 09+O.O6 


287 + !" 5 

-64 


9 9+4.5 
z - -1.1 


8.8/6 


< 0.01 


LOB Seyfert 2s 


Power law 


1 Q 1+0.04 
^ -0.04 






15.9/8 




Cut-off power law 


1 T 02 


JJA+NC 




12.8/7 


0.23 


* Pexrav 


1 Q 1+0.04 




4 +0 5 

u - -0.3 


13.7/7 


0.32 


Pexrav 


1 8+0-2 


^ Z,J -120 


2 +0 / 3 

-0.2 


12.9/6 


0.53 



Notes. NC: the parameter or the 3<r error were not constrained. 
'Pexrav: the high-energy cut off was not included, being fixed to 
its upper limit E c = 10 6 keV. 



Table 5. Results obtained from the spectral analysis of the aver- 
age hard X-ray spectra of PBL and NPBL Seyfert 2s. The incli- 
nation angle i was fixed to 45°. 



Model 


r 


E c 


R 


X 2 /DO¥ 


FTEST 






[keV] 








PBL Seyfert 2s 


Power law 


1 Q4+0.03 
1 ' 7 -0.03 






23.7/8 




Cut-off power law 


1 6 +02 

^ -0.2 






11.7/7 


0.03 


'Pexrav 


1 QS+0.03 
y -0.03 




1 +06 


8.5/7 


<0.01 


Pexrav 


1 9+0.2 


803t 6 c 


7 +a8 

-0.6 


8.5/6 


0.05 


NPBL Seyfert 2s 


Power law 


1 QC+0.09 
1,7 -0.09 






22.9/8 




Cut-off power law 


1 l +05 
-0.5 


49+70 




13.9/7 


0.07 


"Pexrav 


1 Q7+0.09 
1,7 -0.09 




•3+WC 
-2.3 


14.4/7 


0.08 


Pexrav 


1 9 +0 1 


874 +»c 

-548 


9 Q+NC 
z..7_ 23 


13.7/6 


0.21 



Notes. NC: the parameter or the 3cr error were not constrained. 
'Pexrav: the high-energy cut off was not included, being fixed to 
its upper limit E c = 10 6 keV. 



served, one would need a value of the column density as large as 
A^h = 5 x 10 24 cnr 2 . 

We then tested different values of the reflection parameter. 
From the spectral analysis, we found upper limits of the re- 
flection parameter of R < 0.4 for Seyfert Is and Seyfert 1.5s, 
and of R < 0.5 for LOB Sy2s. These values are significantly 
lower than that obtained for MOB Sy2s (R = 2.2^). We sim- 
ulated and compared spectra obtained using Eqs.|3fandS] con- 
sidering different amounts of reflection and absorption for LOB 
(R = 0.5 and R = 0.2, N H = 0, i = 60°) and MOB (R = 1, 
JV H = 7 x 10 23 cnr 2 , 2 = 60° and i = 30°) objects. From Fig.Qj] 
(right panel), one can see that a greater unabsorbed reflection 
component in MOB Seyfert 2s would explain their larger amount 
of flux in the 20-60 keV energy band with respect to LOB ob- 
jects. 

We fitted the spectrum of MOB Sy2s with Eq.|4]in XSPEC, 
setting the column density of the absorbing medium to the av- 
erage value (weighted by the detection significance) of the sam- 
ple (A^h = 3 x 10 23 cm -2 ). We first assumed a edge-on reflector 



(i = 60°), and fixed Ec to different values. The best fit yields 
R > 2 and is obtained for Ec = 300 keV. The value of the re- 
flection parameter is R > 1 even when considering higher values 
of Ec- Considering a face-on reflector (i = 30°), we obtained 
a lower value of the reflection parameter for Ec = 300 keV 
(R = 1.1*H), while we obtained R > 1 for higher values of 
E C - The results of the fits are listed in Table|6] 



7. Discussion 

We now discuss the average hard X-ray spectra of the different 
classes of Seyfert galaxies, the differences and similarities be- 
tween PBL and NPBL Sy2s, the average hard X-ray emission of 
NLSls, and all this in the framework of the unified model and 
the cosmic X-ray background. 



10 



C. Ricci et al.: Reflection in Seyfert galaxies and the unified model of AGN 



Table 6. Results obtained by fitting the spectrum of MOB 
Seyfert 2s with Eq.|4] fixing Nn to the average value of the sam- 
ple. We tested different values of the cutoff energy and the incli- 
nation angle. 



E c 


r 


R 


( 


x 2 


DOF 


[keV] 












300* 


i 7O+0.07 
-0.07 


3 + j° 


60° 


9.6 


7 




i 7 4+0.06 
1 - /H -0.06 


2+3 


45° 


9.7 


7 




1 76 +007 


1 3 +LS 


30° 


9.9 


7 


500* 


1 fi^+0.07 
•° -0.07 


■2 c+yvc 


60° 


10.3 


7 




1 ao+0.06 


2 5 +7 - 8 


45° 


10.3 


7 




1 o 1 +0.07 




30° 


10.4 


7 




1 qq+0.06 
1 - yy -0.06 


1 Q+iVC 

IJ -10 


60° 


11.9 


7 




1 qc+0.06 
7 -0.05 


c o+iVC 

J.J_ 37 


45° 


11.8 


7 




1 96 +007 

i,:7u -0.07 


Q+JVC 

—2 


30° 


11.9 


7 



Notes. /VC: 3<x error not constrained; *: parameter fixed. 



7.1. Seyfert galaxies: the continuum 

From our analysis of the average hard X-ray spectra of Seyfert 
galaxies, we found that the primary power law emission is sim- 
ilar for both Seyfert Is and Seyfert 2s, and can be described by 
a photon index of T ^ 1.8. Comparing the average spectra to a 
simulated power law, we found for Syls, Syl.5s and Sy2s only a 
weak softening in the last energy bin, which might imply that on 
average the cutoff is at energies Ec > 150 keV. This was con- 
firmed by the model-dependent spectral analysis, from which 
we found Ec > 200 keV, which is consistent with the results 
of BeppoSAX observations (Dadina 2008). 

In the thermal Comptonization scenario, the power law 
emission of AGN is commonly assumed to arise from inverse 
Comptonization of soft photons in a hot plasma. Within this sce- 
nario, the cutoff energy and the photon index are linked to the 
temperature and the optical depth of the hot plasma responsi- 
ble for the Comptonization. According to Petrucci et al. (2001), 
the temperature can be estimated as kT e = Ec/2 for r < 1, or 
kT e = Ec/3 for r » 1, and the three parameters are related by 



r - ! 4 9 - + meC ~ 
14 £7yr(i +t/3) 



1/2 



(5) 



Using Eq.[5]and the values of the parameters obtained using the 
pexrav model, we calculated the average optical depths of the 
different classes of Seyferts. For Seyfert Is and Seyfert 1.5s, we 
obtained r = 0.8*q5, which is consistent with the values ob- 
tained for LOB (t = 0.7+° 6) and MQB Sey f ert 2s (r 2 = 0.9^). 
The consistency of r indicates that the average physical charac- 
teristics of the Comptonizing medium are similar for different 
classes of Seyferts, as expected from the UM. 

7.2. Seyfert galaxies: the reflection component 

The hard X-ray spectrum of Seyfert 2s shows a harder emission 
in the 20-60 keV band than those of Seyfert Is and Seyfert 1.5s 
(Fig-0- Most of this difference can be ascribed to the contribu- 
tion of MOB Sy2s (Fig. [9]), while LOB Sy2s show characteristics 
more similar to Seyfert Is and Seyfert 1.5s (Fig. [10}, although 
still having a ~ 20% excess in the 30-50 keV band. 



The greater reflection of MOB Sy2s cannot be explained 
solely by the dampening of the continuum caused by absorp- 
tion, but we have demonstrated that different values of R must 
also be taken into account (Fig. [TBI. We have in fact shown that 
the average value of Seyfert 1/1. 5s and LOB Sy2s is R < 0.4 and 
R < 0.5 (Section|5]l, respectively, while for MOB Sy2s this value 
is greater (R > 1, see Section[6]l. 

From a purely geometrical point of view, and considering the 
disk as the main reflector, one would expect a greater influence 
of the reflection component in the spectrum of Seyfert Is rather 
than in that of Seyfert 2s, the angle i between the normal to the 
disk and the observer being smaller. Different luminosity distri- 
butions could introduce a bias in our hard X-ray selected sample 
of MOB Sy2s, with sources having a stronger reflection compo- 
nent being brighter and thus more easily detected. However, the 
difference between the average luminosities of MOB (L M0B ~ 
1.7 x 10 43 ergs-') and LOB Sy2s (L L0B 2.1 x lO^ergs" 1 ) 
is insignificant. Moreover, a KS test results in a probability of 
86% that the two samples of Seyfert 2s are statistically com- 
patible. Different inclination angles might also play a role in ex- 
plaining the large reflection of MOB Sy2s, although the average 
values of i would have to be extreme to account for the observed 
differences. 

A possible explanation is that in MOB Sy2s the putative 
torus, which is also the dominant reflector, covers a larger 
fraction of the X-ray source than in lightly obscured objects. 
A similar geometry was hypothesized by Ueda et al. (2007) 
to explain the Suzaku spectra of SWIFT J0601. 9-8636 and 
SWIFTJ0138.6-4001. Further evidence for the existence of 
these deeply buried AGN was found by Eguchi et al. (2009). 
Being more absorbed, MOB Sy2s might on average have more 
matter surrounding the active nucleus, which would be responsi- 
ble for the larger reflection. Ramos Almeida et al. (2009, 201 1) 
analyzing the mid-infrared emission of Seyfert galaxies found 
that the clumpy absorbers in Sy2s have larger covering factors 
(C T = 0.95 + 0.02) than those in Syls (C T = 0.5 + 0.1). Among 
the 12 Sy2s in their sample, ten have Nh ^ 10 23 crrT 2 , which 
supports our argument of a larger average covering fraction of 
the absorber for more obscured sources. The small value of the 
average reflection parameter of Syls and Syl.5s might imply 
that on average the reflection from the disk does not play as im- 
portant a role as that of the absorber. 

A value of R > 1 is unphysical if related only to the geom- 
etry, and might imply that part of the direct emission is blocked 
by partially covering material, with the transmission efficiency 
being < l/R (Ueda et al. 2007). If the X-ray source is par- 
tially covered by Compton-thick material, the reflected compo- 
nent would in fact have a stronger relative influence over the 
continuum. Krolik & Begelman (1988) showed that a torus with 
a smooth dust distribution cannot survive close to the AGN, and 
proposed that the material is distributed in a clumpy structure. 
Mid-infrared spectra of Seyfert galaxies have been proven to be 
consistent with the clumpy torus scenario (e.g., Mor et al. 2009). 
Evidence has been found in the past few years of significant vari- 
ations in Nh, with changes from Compton-thick to Compton- 
thin states on timescales from weeks (Risaliti et al. 2005) to 
~ 10 hours (Risaliti et al. 2009). These variations have been in- 
terpreted as eclipses of the X-ray source caused by the BLR. If 
the BLR lies between the X-ray source and the putative torus, 
then it could provide the CT clumps needed to deplete the con- 
tinuum emission. The variations observed in NGC 1365 (Risaliti 
et al. 2005, 2009) might be common in MOB Sy2s, and the ex- 
tremely long INTEGRAL IBIS/ISGRI observations we used to 
derive the average spectra might have registered some of them. 
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The presence of CT clumps might be a common characteristic in 
the absorbers of Seyfert galaxies, and be responsible, at least in 
part, for the large reflection observed in MOB Sy2. 

Both effects might be at work in MOB Sy2s, with more ab- 
sorbed objects having on average more matter around the X-ray 
source, and the distribution of this matter being clumpy. 

7.3. Seyfert 2s with and without PBLs 

Spectropolarimetric surveys indicate that only 30-50% 
Seyfert2s show PBLs (Tran 2001, 2003). The reason for the 
non-detection of PBLs in all Seyfert 2s is still debated, and 
it has been hypothesized (Tran 2001, 2003) that NPBL Sy2s 
might be a different class of Seyfert 2, that possibly lack 
the BLR and have large-scale characteristics more similar to 
Seyfert Is. Deluit (2004) discussed the possible existence of 
differences in the 15-136 keV spectra of Seyfert2s with and 
without PBLs using BeppoSAX data, and concluded that the 
characteristics of PBL Sy2s appear to be more similar to those 
of Seyfert Is than NPBL Sy2s. Using a slightly larger sample, 
but broader energy range and much longer exposures, we 
analyzed the stacked spectra of PBL and NPBL Sy2s, finding 
that the two samples do not show any significant spectral 
differences. The greatest difference between the two samples 
is their luminosities. The average 17-80 keV luminosity of 
PBL Sy2s (L PBL 1.9 x lO^ergs" 1 ) is greater by a factor 
of ^ 3 than that of NPBL Sy2s (L NPBL = 6 x lO^ergs" 1 ). 
The luminosity of the [ OIII] line is often used as an indicator 
of the bolometric luminosity of AGN. Collecting data from 
the literature, we found that the average luminosity of PBL 
Sy2s (L [0 iii] = 4.2 x 10 40 ergs -1 ) is still higher than that of 
NPBL Sy2s (L [0 ni] = lO^ergs" 1 ). Although this does not 
give any indication of the intrinsic luminosity of the BLR, 
it implies that the bolometric luminosity of PBL Sy2s is on 
average greater than that of NPBL Sy2s, confirming what 
was found by Tran (2001) and Lumsden & Alexander (2001). 
We also found that in our sample NPBL Sy2s are on average 
slightly more absorbed (A^ PBL = 3 x 10 23 cnr 2 ) than PBL Sy2s 
(JV PBL = 1.5 x 10 23 cnr 2 ). 

The lack of remarkable differences in the hard X-ray spectra 
of PBL and NPBL Seyfert 2s fits the UM, according to which 
they are supposed to belong to the same family. 

7.4. Narrow-line Seyfert 1s 

The 0.1-10keV spectra of NLSls are steeper than those of their 
broad-line counterparts (e.g., Netzer 2001). Owing to this, few 
of them are detected above 20 ke V, and only a handful of NLS 1 s 
have been studied so far in the hard X-rays. A sample of five ob- 
jects was examined using IBIS/ISGRI data (up to 100 keV) by 
Malizia et al. (2008), who found that their spectra present a con- 
tinuum significantly steeper than Seyfert Is and have an average 
photon index of T = 2.6 + 0.3. They interpreted this as being due 
to a low energy (Ec < 60 keV) of the cutoff. Our sample of 205 
Seyfert galaxies contains 14 NLSls, which makes it the largest 
ever studied at these energies. In our sample, NLSls represent 
15% of the total number of Syls and Syl.5s. This is a much 
smaller fraction than the ~ 50% of ROSAT soft X-ray selected 
NLSls (Grupe 2004), but is consistent with the -15% obtained 
from an optically selected sample (Williams et al. 2002). From 
our analysis, we confirmed that in the hard X-rays the average 
spectrum of NLS 1 s is also steeper than those of Sy 1 s and Sy 1 .5 s, 
with the spectra clearly diverging above ~ 40 keV. This might be 



due to different values of F or to different energies of the cutoff. 
The fact that up to ~ 40keV the average spectrum of NLSls is 
consistent with that of Syls favors the hypothesis that the photon 
indices are consistent, while the energy of the cutoff is lower for 
NLSls. This is supported also by the fact that for the objects for 
which good quality data are available in the soft X-ray band, the 
values of the photon indices are not significantly different from 
those of Syls (Malizia et al. 2008, Jimenez-Bailon et al. 2008, 
Winter et al. 2009). Consistent photon indices might be due to a 
bias towards low values of F for our hard X-ray selected sample. 

7.5. The unified model of AGN 

The fact that the average primary (i.e. excluding the reflected 
component) hard X-ray emission of Seyfert galaxies is indepen- 
dent of their optical classification, and that PBL and NPBL Sy2s 
have consistent hard X-ray spectra fits well the zeroth order UM. 
However, the idea that the anisotropic absorber is the same for 
all classes does not allow us to easily explain the large reflection 
seen in Seyfert 2s, and in particular in MOB Seyfert 2s. 

A possible geometry might be that in MOB Sy2s the absorb- 
ing material covers a large fraction of the X-ray source, which 
is seen through clumps of Compton-thick material, possibly lo- 
cated in the Compton-thin absorber (A^h - 10 23 cm -2 ). This ge- 
ometry would also explain the low hard X-ray to [ OIV] ratio 
found in Compton-thin Seyfert 2s by Rigby et al. (2009), and 
why the luminosity of the Compton-thin Sy2s of our sample is 
significantly lower than that of Syls and Syl.5s. 

7.6. The cosmic X-ray background 

The CXB is known to be produced by the integrated emission of 
unresolved AGN. Thus, to fully understand the spectrum of the 
CXB, it is extremely important to study the average hard X-ray 
spectra of Seyfert galaxies. In their CXB synthesis model, Gilli 
et al. (2007) assumed the disk to be the main reflector, consid- 
ering a weaker reflection for Seyfert 2s than for Seyfert Is, and 
a value of R = 0.37 for the unabsorbed reflector. Using these 
parameters, they hypothesized that, to explain the great amount 
of flux around ~ 30keV, CT Sy2s contribute up to ~ 30% of the 
CXB. Risaliti et al. (1999) and Guainazzi et al. (2005) estimated 
that - 50% of all obscured AGN are CT. However, so far this 
large amount of CT AGN has not been observed (e.g., Ajello et 
al. 2008a, Paltani et al. 2008), and only few of them are known. 
Treister et al. (2009) (hereafter T09) discussed the strong de- 
generacy between the fraction of CT AGN and the value of the 
reflection parameter used, finding that for a value of R ^ 1, a 
lower fraction of CT objects (~ 10%) is needed. Gandhi et al. 
(2007) also showed how an enhanced reflection component, due 
to light bending in their case, might decrease the large fraction 
of heavily absorbed AGN. 

We have shown here that the average spectrum of MOB 
Sy2s has a strong reflection component, and an average value of 
R > 1, although for LOB objects we found a value of R smaller 
than that used by Gilli et al. (2007). This might imply that ob- 
scured objects start to contribute significantly to the peak of the 
CXB already for Nh ~ 10 23 cirT 2 . To test this, we calculated 
the average value of R for AGN, and extrapolated the CT AGN 
density as shown in Fig. 4 of T09. To obtain the average value 
of R, we used the Nh distribution of T09, which considers about 
the same amount of absorbed (AB, Nn > 10 23 crrT 2 , which in- 
cludes MOB and CT AGN) and LOB (N H < 10 23 cnr 2 ) objects. 
We showed that the average spectra of CT and MOB Sy2s are 
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consistent, thus we divided the AGN into two classes, AB and 
LOB, and used their average reflection parameters (Sect. ED to 
calculate the average value of R. The value obtained R ~ 1.1 
agrees with that found by T09, and implies that CT AGN have 
a local density of ~ 2 x 10~ 6 Mpc~ 3 , and represent a fraction of 
~ 10% of the total population of AGN. 

We fitted the INTEGRAL IBIS/ISGRI CXB spectrum of 
Tiirler et al. (2010) with the sum of the best-fit models obtained 
for AB and LOB objectfl fixing all the parameters to the best-fit 
values found in Sect. 15. II leaving only the normalization free to 
vary. Following the Nn distribution of T09, and assuming that all 
the objects contribute equally to the CXB, we forced the two val- 
ues of the normalization to have the same value (Alob = Aab)- 
We set the redshift to z — 1 . We thus obtained a fit with only one 
free parameter, the global normalization. This simplified model 
yields a good fit of the CXB spectrum Or 2 = 13.5 for 15 DOF) 
in the 20-200 keV band. 

In FigJT4j we show our fit to the CXB, and the contribution 
of the different components to the spectrum. Using the A^h distri- 
bution of T09, we transformed the models of the average spectra 
of LOB and AB objects into those of the three classes used by 
Gilli et al. (2007): obscured (10 24 > Na > 10 22 ctrT 2 ), unob- 
scured (A^h < 10 22 cm~ 2 ), and CT objects. We added Compton 
scattering and photoelectric absorption to the model of CT AGN 
(setting A^h = 10 24 cirT 2 ). 



INTEGRAL IBIS/ISGRI 




10 100 
Energy (kc«) 



Fig. 14. Fit to the spectrum of the CXB measured by INTEGRAL 
IBIS/ISGRI (squares, Tiirler et al. 2010). The model used is a 
combination of our best-fit models for the average spectra of the 
different classes of objects (AB and LOB), normalizing their ra- 
tio using the A^h distribution of Treister et al. (2009). We also 
show the data obtained by INTEGRAL JEM-X (stars, Tiirler et 
al. 2010) and Swift BAT (triangles, Ajello et al. 2008b) above 
lOkeV. 



8. Summary and conclusion 

We have presented a study of the average hard X-ray spectra of 
Seyfert galaxies performed using INTEGRAL IBIS/ISGRI data. 
In the following, we summarize our findings: 

• Different classes of Seyfert galaxies have the same average 
nuclear continuum (T 1.8, Eq > 200 keV), and their 
Comptonizing medium has on average a consistent value of 
the optical depth (t ^ 0.8). 



• NLSls have a steeper hard X-ray spectrum than Seyfert Is 
and Seyfert 1.5s above ~ 40 keV. This is likely due to a lower 
cutoff energy, as proposed by Malizia et al. (2008). 

• PBL and NPBL Sy2s have similar hard X-ray spectra, the 
only difference between the two classes in our incomplete 
sample being their luminosity. 

• Seyfert 2s have on average a stronger reflection component 
than Seyfert Is and Seyfert 1.5s. 

• Most of the reflection of Seyfert 2s comes from MOB 
(10 23 cnr 2 < N H < 10 24 cnr 2 ) Sy2s. 

• The large amount of reflection observed in MOB Seyfert 2s 
might be explained by an X-ray source highly covered by the 
absorbing material. The Compton-thin absorber might con- 
tain CT clumps, partially covering the X-ray source. 

• The large amount of reflection of MOB Sy2s reduces by a 
factor of ~ 3 the amount of CT Seyfert 2 needed to explain 
the CXB peak. Our results are consistent with the fraction of 
CT AGN being - 10%. 

Our results, while confirming that the X-ray engine is the same 
for all classes of Seyfert galaxies, point towards the existence of 
significant differences in the structure of the medium surround- 
ing the X-ray source between LOB Seyfert galaxies and MOB 
Seyfert 2s. These differences cannot be solely due to viewing 
angle dependence, but might imply that there are morphologi- 
cal differences between LOB and MOB AGN. 

Appendix A: List of sources 

In the following, we list the coordinates, redshifts, detection sig- 
nificances, count rates, column densities, and luminosities of the 
samples of Seyfert Is (Table lATTT i. Seyfert 1.5s (1A.2I >, Compton- 
thin Seyfert 2s COi CT Seyfert 2s SKM . and NLSls (|A"3T > used 
for the stacking analysis. 

Appendix B: Spectra 

In Fig. IB. II we show the average 17-250keV spectra of 
Seyfert Is, Seyfert 1.5s, Seyfert2s, CT Seyfert2s, and NLSls. 
The parameters of the models used for the fit are reported in 
Table|3] In Fig |B.2l we show the average spectra of MOB and 
LOB Seyfert 2s. The parameters of the models are reported in 
Table@ 
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Table A.l. Coordinates, redshifts, detection significances, count rates, column densities, and luminosities of the Seyfert Is used for 
the stacking analysis. Detection significances, count rates, and luminosities are in the 17-80 keV band. The flux of the Crab in the 
17-80 keV band is of 278 ct s- 1 . 



Name 


RA 


Dec 


z 


Det. significance 


Count rate 


N H 


logLl7-80keV 




[h m s] 


[° ' "] 




M 


r * -In 
[cts 'J 


[10 cm z \ 


r -In 


IGR J02086-1742 


02 08 31.4 


-17 43 04.8 


0.129 


7.3 


0.31 ±0.04 


< 0.017° 


44.78 


IGR J02097+5222 


02 09 37.7 


+52 26 43.6 


0.049 


13.5 


0.56 + 0.04 


0.03 6 


44.20 


Mrk 590 


02 14 33.6 


-00 46 00.3 


0.027 


5.0 


0.13 ±0.03 


0.03 c 


43.05 


Mrk 1040 


02 28 14.6 


+31 18 39.4 


0.016 


8.4 


0.72 ± 0.09 


0.067 rf 


43.33 


RBS 345 


02 42 16.0 


+05 31 48.0 


0.069 


5.1 


0.23± 0.04 


_ 


44.11 


ESO 548-01 


03 41 54.7 


-21 15 28.8 


0.014 


6.1 


0.64 ±0.11 


_ 


43.17 


UGC 3142 


04 43 46.9 


+28 58 19.0 


0.022 


13.8 


0.58 ± 0.04 


1.4" 


43.52 


LEDA 168563 


04 52 04.7 


+49 32 45.0 


0.029 


11.6 


0.51 ±0.04 


< 0.22* 


43.70 


Ark 120 


05 16 11.5 


-00 09 00.6 


0.033 


29.3 


0.77 ± 0.03 


< 0.1 c 


43.99 


EXO 055620-3820.2 


05 58 02.0 


-38 20 01.0 


0.034 


5.6 


0.31 ±0.06 


2.57 / 


43.62 


IRAS 05589+2828 


06 02 09.7 


+28 28 17.0 


0.033 


17.5 


0.61 ±0.03 


< 0.04 6 


43.89 


SWIFT J0640.4-2554 


06 40 10.8 


-25 49 48.0 


0.026 


5.0 


0.39 ± 0.08 


0.33^ 


43.49 


IGR J07597-3842 


07 59 41.7 


-38 43 57.4 


0.040 


21.8 


0.60 ± 0.03 


0.05 rf 


44.05 


PG 0804+761 


08 10 58.7 


+76 02 42.5 


0.100 


5.0 


0.20 ± 0.04 


0.023 rf 


44.37 


Fairall 1146 


08 38 30.7 


-35 59 35.0 


0.032 


14.0 


0.34 ± 0.02 


0.l d 


43.61 


IGR J09026-4812 


09 02 37.3 


-48 13 34.1 


0.039 


21.6 


0.37 ± 0.02 


0.9* 


43.82 


SWIFT J0917.2-6221 


09 16 09.4 


-62 19 29.5 


0.057 


12.1 


0.34 ± 0.03 


0.5 h 


44.11 


4U 0937-12 


09 45 42.1 


-14 19 35.0 


0.008 


17.1 


0.99 ± 0.06 


1.19 f 


42.87 


SWIFT J1038.8-4942 


10 38 45.0 


-49 46 55.0 


0.060 


9.3 


0.25 ± 0.03 


1.55 / 


44.02 


IGR Jl 1457- 1827 


11 45 41.0 


-18 27 29.0 


0.033 


6.8 


0.87 ±0.13 


^ 


44.05 


IGR J12107+3822 


12 10 43.4 


+38 22 51.6 


0.023 


5.3 


0.17 ±0.03 




43.02 


IGR J12136-0527 


12 13 37.0 


-05 26 53.0 


0.066 


5.3 


0.16 ±0.03 


2.0' 


43.91 


IGR J12172+0710 


12 17 09.0 


+07 09 33.0 


0.007 


8.3 


0.22 ± 0.03 


0.15* 


42.10 


Mrk 50 


12 23 24.1 


+02 40 44.8 


0.023 


8.9 


0.22 ± 0.02 


0.018 rf 


43.14 


NGC 4593 


12 39 39.4 


-05 20 39.3 


0.009 


39.7 


1.03 ±0.03 


0.02^ 


42.99 


ESO 323-77 


13 06 26.6 


-40 24 50.0 


0.015 


16.6 


0.50 ± 0.03 


6.0* 


43.12 


IGR J13109-5552 


13 10 43.1 


-55 52 11.7 


0.104 


15.2 


0.31 ± 0.02 


< o.r 


44.60 


Mrk 279 


13 53 03.5 


+69 18 29.2 


0.031 


5.0 


0.64 ±0.13 


0.02^ 


43.86 


RHS 39 


14 19 22.2 


-26 38 41.0 


0.022 


10.4 


0.46 ± 0.04 


< 0.05 r 


43.42 


IGR J14471-6414 


14 46 28.3 


-64 16 24.1 


0.053 


9.7 


0.19 ± 0.02 


< O.V 


43.80 


IGR J161 19-6036 


16 11 51.4 


-60 37 53.1 


0.016 


14.8 


0.33 ± 0.02 


0A d 


43.00 


IGR J16482-3036 


16 48 14.9 


-30 35 06.1 


0.031 


29.9 


0.62 ± 0.02 


0.13 c 


43.85 


IGR J16558-5203 


16 56 05.7 


-52 03 41.2 


0.054 


23.9 


0.47 ± 0.02 


o.oir 


44.21 


TGR 11741 8-1212 


17 41 55 3 


-12 11 57.5 


0.037 


19.4 


43 + 02 


0.1 rf 


43.84 


IGR J18027- 1455 


18 02 48.0 


-14 54 54.8 


0.035 


23.6 


0.42 ± 0.02 


0.41' 


43.78 


IGR J 18259-0706 


18 25 57.6 


-07 10 22.8 


0.037 


13.2 


0.24 ± 0.02 


0.6 C 


43.58 


ESO 140-43 


18 44 54.0 


-62 21 52.9 


0.014 


11.8 


0.52 ± 0.04 


1.8" 


43.08 


IGR J18559+ 1535 


18 56 00.6 


+ 15 37 58.0 


0.084 


18.3 


0.35 ± 0.02 


0.7" 


44.46 


ESO 141-55 


19 21 14.2 


-58 40 15.0 


0.037 


12.7 


0.64 ± 0.05 


0.004 r 


44.01 


SWIFT J1933.9+3258 


19 33 47.3 


+32 54 25.0 


0.056 


7.0 


0.22 ± 0.03 


< 0.04 m 


43.91 


IGRJ19405-3016 


19 40 15.2 


-30 15 48.5 


0.052 


9.0 


0.29 ± 0.03 


< 0.1 c 


43.97 


RX J2135.9+4728 


21 35 54.4 


+47 28 28.3 


0.025 


14.1 


0.28 ± 0.02 


0.4 C 


43.31 


IGR J22292+6647 


22 29 13.5 


+66 46 51.8 


0.113 


9.7 


0.22 ± 0.02 


0.2 C 


44.52 


IGR J23 206+6431 


23 20 36.8 


+64 30 42.8 


0.072 


9.0 


0.15 ±0.02 


0.6 C 


43.96 



Notes. " Rodriguez et al. (2010); b Winter et al. (2008) and references therein; c Beckmann et al. (2009) and ref. therein; 
d Bodaghee et al. (2007) and ref. therein; e Ricci et al. (2010); f Winter et al. (2009); s Tomsick et al. (2008); h Malizia 
et al. (2007); ' Paltani et al. (2008) and ref. therein; ; Beckmann et al. (2006) and ref. therein; * Jimenez-Bailon 
et al. (2008); ' Panessa et al. (2008); m Landi et al. (2007a). 
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Table A.2. Coordinates, redshifts, detection significances, count rates, column densities, and luminosities of the Seyfert 1.5s used 
for the stacking analysis. Detection significances, count rates, and luminosities are in the 17-80 keV band. The flux of the Crab in 
the 17-80 keV band is of 278 ct s- 1 . 
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Table A.3. Coordinates, redshifts, detection significances, count rates, column densities, and luminosities of the Seyfert 2s used for 
the stacking analysis. Detection significances, count rates, and luminosities are in the 17-80 keV band. The flux of the Crab in the 
17-80 keV band is of 278 ct s- 1 . 
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Notes. References for the column density. " Landi et al. (2007b); * Beckmann et al. (2009) and ref. therein; ' Bodaghee et al. (2007) 
and ref. therein; d Pfeerkorn et al. (2001); ' Masetti et al. (2008); 1 Rodriguez et al. (2010); * Greenhill et al. (2008); * Winter et al. (2009) 
and ref. therein; ' Landi et al. (2010); > Winter et al. (2008) and ref. therein; ' Paltani et al. (2008); ' Risaliti et al. (2000); 
"' Rodriguez et al. (2008); " Dadina et al. (2007); " Bikmaev et al. (2008). 

References for the PBL. 1 Shu et al. 2007 and ref. therein; 2 Veron-Cetty & Veron 2010 and ref. therein; 3 Cai et al. 2010; 
4 Nagao et al. 2000; 5 Barth et al. 1999. 
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Table A.4. Coordinates, redshifts, detection significances, count rates, column densities, and luminosities of the Compton-thick 
Seyfert 2s used for the stacking analysis. Detection significances, count rates, and luminosities are in the 17-80 keV band. The flux 
of the Crab in the 17-80 keV band is of 278 cts" 1 . 



Name 


RA 


Dec 


z 


Det. significance 


Count rate 




logLl7-80keV 




[h m s] 


[° ' "] 




M 


r . -In 

[cts '] 


[10 z/ cm z ] 


r -In 

Lergs '] 


NGC 1068 


02 42 40.8 


-00 00 48.4 


0.003 


14.0 


0.45 ± 0.03 


150" 


41.91 


SWIFT J0601.9-8636 


06 05 39.1 


-86 37 52.0 


0.006 


7.3 


0.40 ± 0.06 


101* 


42.40 


Mrk3 


06 15 36.3 


+71 02 14.9 


0.014 


31.3 


1.43 + 0.05 


110 c 


43.70 


NGC 3281 


10 31 52.1 


-34 51 13.3 


0.011 


12.0 


0.78 ± 0.07 


I5l d 


43.26 


IGRJ13091+1137 


13 09 05.6 


+ 11 38 02.9 


0.025 


14.8 


0.51 ± 0.03 


90 c 


43.73 


NGC 5643 


14 32 41.3 


-44 10 24.0 


0.003 


9.2 


0.20 ± 0.02 


60 - 100 e 


41.50 


Mrk 477 


14 40 38.1 


+53 30 15.9 


0.037 


5.8 


0.44 ± 0.07 


> ioof 


44.09 


NGC 5728 


14 42 23.3 


-17 16 51.6 


0.009 


13.6 


0.80 ± 0.06 


210° 


43.12 


IGRJ16351-5806 


16 35 13.2 


-58 04 49.7 


0.009 


12.9 


0.28 ± 0.02 


> 150 g 


42.69 


NGC 7582 


23 18 23.5 


-42 22 14.1 


0.005 


5.5 


1.20 + 0.20 


160" 


42.77 



Notes. 5 Delia Ceca et al. (2008) and ref. therein; * Ueda et al. (2007); c Beckmann et al. (2006) and ref. therein; 
d Bodaghee et al. (2007) and ref. therein; e Guainazzi et al. (2004); f Bassani et al. (1999); * Malizia et al. (2009). 



Table A.5. Coordinates, redshifts, detection significances, count rates, column densities, and luminosities of the Narrow Lines 
Seyfert Is used for the stacking analysis. Detection significances, count rates, and luminosities are in the 17-80 keV band. The flux 
of the Crab in the 17-80 keV band is of 278 ct s" 1 . 



Name 


RA 


Dec 


z 


Det. significance 


Count rate 




login-SOkeV 




[h m s] 


[° ' "] 




M 


[cts- 1 ] 


[10 22 cm" 2 ] 


[ergs -1 ] 


1H 0323+342 


03 24 41.2 


+34 10 45.9 


0.063 


5.9 


0.20 ± 0.03 


0.1" 


43.97 


Mrk 1 10 


09 25 12.9 


+52 17 10.5 


0.035 


5.0 


0.86 + 0.17 


0.019° 


44.09 


NGC 4051 


12 03 09.6 


+44 31 53.2 


0.002 


13.6 


0.47 ± 0.03 


< 0.02* 


41.34 


Mrk 766 


12 18 26.6 


+29 48 45.6 


0.013 


7.7 


0.26 ± 0.03 


0.8 C 


42.71 


NGC 4748 


12 52 12.0 


-13 25 48.0 


0.014 


5.2 


0.16 + 0.03 


-0 d 


42.57 


Mrk 783 


13 02 58.8 


+ 16 24 27.5 


0.067 


7.3 


0.25 ± 0.03 


0.046° 


44.12 


NGC 5506 


14 13 14.9 


-03 12 27.0 


0.007 


47.2 


2.7 ± 0.06 


3.4* 


43.19 


IGRJ14552-5133 


14 55 17.9 


-51 34 13.3 


0.016 


12.5 


0.25 ± 0.02 


0.1° 


42.88 


IRAS 15091-2107 


15 11 54.0 


-21 21 21.6 


0.044 


5.3 


0.22 ± 0.04 


0.7 P 


43.70 


IGRJ16185-5928 


16 18 26.4 


-59 28 45.3 


0.035 


10.7 


0.23 ± 0.02 


< 0.P 


43.52 


IGRJ16385-2057 


16 38 31.1 


-20 55 25.0 


0.027 


7.7 


0.26 ± 0.03 


0.21* 


43.35 


IGRJ19378-0617 


19 37 33.1 


-06 13 04.0 


0.010 


11.2 


0.32 ± 0.03 


=a 0* 


42.58 


ESO 399-20 


20 06 57.2 


-34 32 54.0 


0.025 


6.5 


0.23 ± 0.04 


0.048° 


43.23 


IGR J21277+5656 


21 27 45.9 


+56 56 34.4 


0.014 


24.2 


0.52 ± 0.02 


0.1° 


43.08 



Notes." Bodaghee et al. (2007) and ref. therein; Beckmann et al. (2009) and ref. therein; 

Beckmann et al. (2006) and ref. therein; d Grupe et al. (1998); e Jimenez-Bailon et al. (2008); 

f Malizia et al. (2008); « Rodriguez et al. (2008),; * Malizia et al. (2008). 
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Fig. B.l. Spectra of different classes of AGN obtained using the pexrav model. We also show the contributions to the^ 2 relative to 
the best fits. The results of the spectral analysis are reported in Table|3] 



C. Ricci et al.: Reflection in Seyfert galaxies and the unified model of AGN 

MOB Seyfert 2 LOB Seyfert 2 



19 



20 



50 



100 



200 



Energy (keV) 



20 



50 100 
Energy (keV) 



200 



Fig. B.2. Spectra of MOB {left panel) and LOB {right panel) Seyfert 2s, the model used is a cutoff power law plus a reflection 
component from neutral matter ipexrav). The results of the analysis are reported in Tabled 



